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Abstract. 

We present the first completely analytical computation of 
the full differential 7-7 pair production rate from compact ra- 
diation fields, exact to 2nd order QED, and use this result to 
investigate the validity of previously known approximations. 

Key words: plasmas — radiation mechanisms: Pair produc- 
tion — gamma-rays: theory 



1. Introduction 

The discovery of high-energy 7-radiation from extragalactic 
compact objects has motivated many authors to consider the 
effects of 7-ray absorption by 7-7 pair production, eventually 
inducing pair cascades. The relevance of 7-7 pair production 
to astrophysical systems has first been pointed out by Nikishov 
(1962). The first investigation of the 7-7 absorption probability 
of high-energy photons by different soft photon fields, along 
with some useful approximations, can be found in Gould & 
Schreder (1967). 

The energy spectrum of injected electrons and positrons 
due to this process has been studied by several authors (e. g., 
Bonometto & Rees 1971, Aharonian et al. 1983, Zdziarski & 
Lightman 1985, Coppi & Blandford 1990). In most astrophysi- 
cally relevant cases, simple approximations can be used for this 
purpose, without much loss of accuracy. These usually rely on 
the high-energy photon having much higher energy than the 
soft photons and thus dominating the energy input and de- 
termining the direction of motion of the center-of-momentum 
frame of the produced pairs. Bonometto & Rees (1971) used 
basically the same technique as we do, but restricted their anal- 
ysis to the case £1 3> £2, and did not solve the problem analyt- 
ically. Two recipes to calculate the full energy-dependence of 
the injected pairs have been published (Aharonian et al. 1983 
and Coppi & Blandford 1990), but here the reader is still left 
with integrations to be carried out numerically. 

It is the purpose of this paper to derive the full energy- 
spectrum of pairs, injected by 7-7 pair production, exact to 
second order QED for the case of isotropic radiation fields. In 
Section 2, we give a short overview of the kinematics which 
are used in Section 3 to calculate the pair injection spectrum. 



In Section 4, we compare our results to well-known approxi- 
mations and specify the limitations of the various approxima- 
tions. Our analysis is easily generalized to non-isotropic radi- 
ation fields. The derivation presented here is widely analogous 
to the derivation of the pair annihilation spectrum, given by 
Svensson (1982). 

2. Kinematics 

We consider an isotropic photon field n p h(e) where e = 
hi> / (m e c 2 ) is the dimensionless photon energy in a rest frame 
which we call the laboratory frame. The Lorentz invariant 
scalar product of the four-momenta £1,2 of two photons having 
energies £1,2 colliding under an angle of cosine n — cos0 M in 
the laboratory frame is then given by 



£1 • £2 = £162(1 — At) = 2£ c 



(1) 



Here, £ cm is the photon energy in the center-of-momentum 
frame. In order to allow for the possibility to create an electron- 
positron pair, conservation of energy implies e cm = 7 cm , and 
the condition 7 cm > 1 determines the pair-production thresh- 
old. 7 cm is the Lorentz factor of the electron/positron in the cm 
frame where the produced electrons move with speed ±/3 cm c 
and /3cm = yl- 1 /7? m • The definition of the angle variables 
needed in this calculation is illustrated in Fig. 1. 
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Fig. 1. Definition of the angles in cm and laboratory frame. ei : 2 
denotes the direction of motion of an incoming photon, f3 C m,± is the 
direction of motion of the produced electron and positron in the cm 
and the laboratory frame, and /3 C characterizes relative motion of 
the laboratory and the cm-frame, respectively. 
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The cm frame moves relative to the laboratory frame with 
velocity — (3 c c and Lorentz factor f c = (1 — /3 2 ) -1 / 2 . The four 
velocity of the laboratory frame (/3 c c in the cm frame) is de- 
noted by vl_. The Lorentz factors of the produced pairs in the 
laboratory frame are related to the cm quantities by 

7± = 7cm7c (1 ± I3 cm ,l3 c u). (2) 

Evaluating the Lorentz invariant scalar product 
€1,2 • VL_ = €1,2 = tcm 7c(l ± PcZ) (3) 
in the laboratory and the cm-frame, respectively, we find 
E 



7c 



2 tcm 

and 



where E — ei + £2, 



(4) 



z — cos9 z = —Jj—y where N = \/ E 2 — 4 eL- (5) 

Inserting Eq. (4) into Eq. (2) and using energy conservation 
{(-cm = 7cm) fixes the angle cosine u to 



M = MO 



E-2j- 

PcmN 



(6) 



The differential cross section for 7-7 pair production (see 
Eq. [11]) depends on 



x = cos6 x — uz + \J\ — u 2 \J\ — z 2 cos 4>- 
3. The pair yield 

The differential yield of produced pairs is calculated as 

Mi-) 



(7) 



00 00 ei* 



da 



dem p h{ei) j de2n ph (e2) / rf^ (1 — n) (8) 



max| Jj- , 7_ +1 — ei [ 



where 



da _ /" ^ d 2 a d-y cm 

dj- J ° m dQcmdjcm dj- 



(9) 



The differential cross section has been evaluated by Jauch 
& Rohrlich (1959): 
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dQcmd'y 

err 

where 



^(^cm 7cm) 
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dn c rn 



1 3 f3c 



(10) 



-1 + 



3 - I3t 



1 (3 cm X 1 H - fjcmX 



(11) 



We may express the solid angle element dQ cm = dudcj>. 
Using Eq. (2), we find 



e/ \ d^fcm j-/ \ ^ 

<K e cm — 7cm) —j = 0(U — Wo) ■ 



C?7_ 



(12) 



This enables us to carry out the w-integration in Eq. (9) 
immediately. If we write the denominators in Eq. (11) as 

1 ± ficmX — a± +b± cos <f> (13) 

with 

a± = l±uoz/3 cm , b± = ±y/l - u 2 Q \Jl - z 2 l3 cm , (14) 
we find 

da 3 1 



277 



3 - a* 

1 -\- ' 



dj- 8 Ne; 

1 



2 

cm 



+ 



2 V a + + b + cos <f> a- + b- cos 



1 1 

+ 



2 eim \ [a+ + b + cos <p] 2 [a- + b- cos 4>] 2 
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= -a T — 

O tcTT) 



1 , 3-/? c 4 



+ 



iVtcm 4 

where 



■ (G+ + G- 



8e? 



G± 



V / eie 2 _ + 
rf± - 2 e; 



2 

^cm C± 



2 

cm 



(eie 2 + elmC±) 3/2 ' 
with 

c± = (ei, 2 - 7-) 2 - 1, 
d± = £1,2 + ei£2 ± 7- ( £ 2 — ei) 
and we used the integrals 



2tt 



a + bcos0 V«2 _ b 2 



d<f> 



2n a 



(a + bcosfi) 2 (a 2 ~6 2 ) 3 / 2 



and the identity 



2 ,2 /i£2 + e cm c± 

a ± - 6 ± = 4 — 



(F + + F_)|> (15) 

(16) 
(17) 

(18) 
(19) 

(20) 
(21) 



(22) 
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which follows from Eqs. (5), (6) and (14). Now, inserting 
Eq. (15) into Eq. (8) yields the exact expression for the dif- 
ferential pair injection rate. Using Eq. (1) we transform the /i 
integration into an integration over de cm . This leads us to 



For c± — we find 



d± 



8 J (eie 2 ) 3 / 2 



= 7 a r c / dei n ph {e-C) 



4e c 



(29) 



OO 

/ 



dt2 n p h{t2) 



\ Ne, 



max { 7f ' T-+ 1_e l } 

3 -A 4 



(G+ + G-)- 



(F+ + F.)\ (23) 



which can be calculated analytically. The integration limits 
follow from 1 < e cm < v /eie2 and the condition 7_ <7 cm 7 c (l± 
PcmPc) which yields 




- ^ I5 £cm } 



(24) 



\ (l-[E - 7-] + 1 ± ^(j-[E-"/-] + l) 2 -E^ 



(25) 



Using the integrals 2.271.4, 2.271.5, 2.272.3, 2.272.4, and 
2.275.9, of Gradshteyn & Ryzhik (1980), wc find as final result 
for the differential pair yield 



. , . 3 f n p h(ei) 
nil-) = 4 a Tcj de 1 — — 



d€2 



n vh (t2) 



<{-L, 7 _+i- ei } 



(26) 



where for c± / we have 
H± 



d± + 2_ 



+1 2 -^"* 



V e i £ 2 + c±e; 



2 / f I 

cm I £ cm ^ 



c± £ cm £i£2 



(27) 



and 



1 



|= In [e cm ^/c± + \Jei£2 + c±e 2 cm j for c± > 
arcsin ^m^-^j for c± < 0. 



(28) 



4. Comparison to approximations 

Now, we use the exact expression, given in Eq. (26) to specify 
the regimes of validity and the limitations of various approxi- 
mations. The first detailed computation of the pair production 
spectrum was presented by Bonometto & Rees (1971). Based 
on the neglect of the energy input of the soft photon, they ba- 
sically follow the same procedure as described above, but do 
not carry out the angle-integration (integration over e cm in our 
formalism) analytically. In the case €1 ^ £2, it is in very good 
agreement with the exact result, but its evaluation is even more 
time-consuming than using the latter. For this reason, we will 
not consider it in detail, but concentrate on approximations 
which really yield simpler expressions than the exact one. 

4-1. S -function approximation for power-law spectra 

Probably the simplest expression for the pair spectrum injected 
by 7-rays interacting with a power-law of energy spectral index 
a' = a — 1 (71(62) oc ^2") i s based on the assumption €1 >• £2 
and on the well-known fact that photons of energy ei inter- 
act most efficiently with photons of energy £2 = 2/ei which 
motivates the approximation for the 7-7 opacity of Gould & 
Schroder (1967) using a S function approximation for the cross 
section, 

ir 77 (ei,e2) ~ -o- T e28 ^e 2 — — ^ 

Since in this approach pair production takes place only near 
the pair production threshold (7 cm w 1), the produced pairs 
have energies 7+^7-^ (ei+e2)/2 (Bonometto & Rees 1971). 
The resulting pair injection spectrum is therefore 



■ 1 \ ~ 1 i\ n p h(eo) 
n(7) « n{a ) ca T n ph 



(5) 



(30) 



(e. g. Lightman & Zdziarski 1987) where eo = 7+ -\/7 2 — 1 
(which, of course, reduces to 27 for 7 3> 1) and r}(a') is a nu- 
merical factor, depending only on the spectral index of the soft 
photon distribution. This approximation yields useful results, 
if the power-law photon spectra extend over a sufficiently wide 
range (eu/eL >,10 4 ) and if for every high-energy photon of en- 
ergy ei there is a soft photon of energy £2 = 2/ei. Else, the in- 
jection spectrum calculated with Eq. (30) cuts off at the inverse 
of the respective cutoff of the soft photon spectrum, seriously 
underpredicting the injection of pairs of higher or lower energy, 
respectively, where the injection spectrum declines smoothly. 
Nevertheless, these pairs can still carry a significant fraction 
of the injected power. Eq. (30) fails also to describe the in- 
jection of low-energetic pairs in case of a high lower cut-off of 
the 7-ray spectrum even if soft photons of energy £2 = 1/ei 
are present. For example, in the case of the interaction of a 
power-law spectrum extending from £1 = 10 2 - 10 6 with a soft 
power-law spectrum extending from £2 = 10~ 7 - 10~ 2 Eq. (30) 
overpredicts the injection of pairs slightly above 7 = 50 by an 
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order of magnitude and cuts off below this energy. A similar 
problem arises at the high-energy end of the injection spec- 
trum. An example for this fact is shown in Fig. 2. In contrast, 
the approximation (30) can well be used to describe the in- 
jection of pairs at all energies if both photon fields extend up 
to (and down to, respectively) £1,2 ~ 1. For more general soft 
photon distributions which are different from a power-law (e. g. 
a thermal spectrum) the analogous 5-function approximation 
has first been introduced by Kazanas (1984). 
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Fig. 2. Differential pair injection rate (arbitrary units) for the in- 
teraction of a power-law from ei = 10 2 - 10 6 , photon spectral index 
ct 7 = 1.5, with a soft power-law from £2 = 10~ 7 - 10~ 2 , a s = 1.5. 



4-2. S -function in electron energy 

Using the full cross section for 7-7 pair production as given 
by Jauch & Rohrlich (1959) instead of the <5-function approx- 
imation adopted in Eq. (30) does not reduce the limitations 
of this power-law approach significantly, but other soft pho- 
ton distributions can be treated more successfully with this 
approximation which in the limit £1 3> £2 reads 



^(7) w 2 c 



ph 



(27) 



(27) 2 



de-2 



(2) 
ph 



(£2) 



(ft) 



ds S<7 77 (s) 



(31) 



1/27 



where s 



and the limits e 



U/L 



are given by Eqs. (24) 
and (25). Here, we have assumed that the produced electron 
and positron have energy 7- = 7+ = £i/2. This approach 
works equally well for power-law photon fields, but in con- 
trast to Eq. (30), it tends to underpredict the injection of low- 
energetic pairs. The same is true for the interaction of 7-ray 
photon fields with thermal soft photon fields where the high- 
energy tail of the injection spectrum is described very accu- 
rately (a few % error) by Eq. (31). The accuracy of this ap- 
proximation improves with decreasing lower cut-off of the 7-ray 
spectrum. E. g., the injection due to a power-law 7-ray spec- 
trum from ei = 2- 10 6 interacting with a thermal spectrum of 
normalized termpature O = 0.1 is described by Eq. (31) with 
a deviation of less than 30 % from the exact result down to 
7 ~ 1. For 7 > 2, the deviation was much less than 10 %. We 
show an example for the latter situation in Fig. 3. 
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Fig. 3. Differential pair injection rate (arbitrary units) for the inter- 
action of a power-law 7-ray spectrum from e\ 
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5 to 10 6 , a-y = 1.5, 
with a thermal blackbody spectrum of temperature — m -2 



4-3. Approximation by Aharonian et al. 

A very useful approximation to the pair injection spectrum 
for all shapes of the soft photon spectrum under the condi- 
tion £2 <C l<ei has been found by Aharonian et al. (1983). 
They use a different representation of the pair production 
cross section and end up with a one-dimensional integral over 
k = yef + £j + 2£i£2/i which is equivalent to our e cm integra- 
tion in Eq. (21). They solve this integration analytically after 
simplifying the integrand and the integration limits according 
to the assumptions mentioned above. The resulting injection 
spectrum is 



71(7) 



32 



CUT 



C?£2 



*7(«1 — f) 



I 4£ 2 ln (^^-^)_ 8eie2 



+ 



2(2£i£ 2 - l)e{ 
7(ei - 7) 



(*--) 

V £l£2/ 



7 2 (£i - 7) 2 



(32) 



It describes the power-law tail of the pair spectrum in- 
jected by power-law 7-ray photon fields perfectly and is much 
more accurate to the injection of low-energy pairs. Interaction 
with a power-law soft photon field is reproduced within errors 
of only a few %. Even if as well the 7-ray as the soft photon 
spectrum extend to £1,2 = 1, the error at 7 ~ 1 increases only 
to ~ 10 %. Problems with this approximation arise if the soft 
photon spectrum extends up to £2 ~ 1, but the 7-ray spec- 
trum has a lower cut-off £1 >• 1. In this case, the injection 
of low-energetic pairs is seriously overpredicted by Eq. (32). 
For power-law soft photon fields, the integration over £2 in Eq. 
(32) can be carried out analytically, as was found by Svensson 
(1987). His Equation (B8) (multiplying with the total absorp- 
tion coefficient oc e a by which the total injection rate had been 
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normalized to 1) is in perfect agreement with the numerical 
results according to Eq. (32). 

The interaction of power-law 7-ray spectra with thermal 
soft photon fields is generally described within an error of a 
few % at all electron/positron energies if the soft photon tem- 
perature is 0<O.l, even if the 7-ray spectrum extends down to 
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Fig. 4. Differential pair injection rate (arbitrary units) for the in- 
teraction of a thermal blackbody spectrum of temperature © = 5 
with itself. 



Interestingly, even the interaction of a mildly relativistic 
thermal photon field (0^<3) with itself (for which Aharonian's 
approximation was not designed) is reproduced reasonably 
well, but the result of Eq. (32) differs from the exact injection 
rate by a roughly constant factor. When artificially introducing 
a factor adjusting the high-energy tails of the injection spectra, 
Eq. (32) overpredicts the injection of low-energetic pairs by a 
factor of ~ 3, but for 7^3 there is very good agreement with 
the exact result. The deviation becomes more important with 
increasing photon temperature, and for = 5 the injection of 
cold pairs is already overpredicted by a factor of > 10. Fig. 
4 illustrates the accuracy of the various approximations for a 
compact thermal radiation of temperature 6 = 5. 

We find that all the statements on soft photon or 7-ray 
power-law spectra made above are only very weakly dependent 
on the respective spectral index. 



5. Summary and conclusions 

We derived the full energy spectrum of injected pairs, produced 
by 7-7 pair production and compared the result to the vari- 
ous approximations known before. We found that the simplest 
expressions, based on S function approximation to the cross 
section can well reproduce the power-law tail of the injection 
spectrum, but have problems at low particle energies. 

The approximation of Aharonian et al. (1983) was shown 
to be the most accurate one and even yields useful results in 
regimes for which this approach was not designed. 
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